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Abstract 

We calculate the strong radiative corrections to the process pp — > W(-^ ev)g*(-^> 
bb). At the Tevatron this process is the largest background to the associated Biggs 
boson production pp — > W(— > ev)H{-^ bb). The calculation is based on the sub- 
traction procedure, and the corrections are found to be significant. 



1 Introduction 



The search for the origin of electroweak symmetry breaking is at the top of the agenda for 
the next generation of collider experiments. In the standard model (SM) the mechanism of 
electroweak symmetry breaking predicts the existence of a single uncharged Higgs boson 
whose mass is a priori unknown. It is currently constrained to be 90 < Mh < 280 GeV 
at 95% confidence level where the upper bound comes from precision electroweak 
measurements, while the lower bound is determined by direct searches.^] 

At the Tevatron collider there is the potential to look for the SM Higgs boson using 
the decay mode H —>■ bb [|J. In the lowest order (LO) the most promising process is 

pp -> W(-> ev)H{-^ bb) . (1) 

This search can cover the mass range up to about 130 GeV, once event samples, perhaps 
as large 30fb _1 , have been accumulated. Collection of data samples of this size will not be 
easy. However, the Tevatron search is of great importance, especially because the mass 
range between 100 < Mh < 130 GeV is one of the most challenging regions for the LHC 
to look for the SM Higgs [§. 

In this letter we calculate the strong radiative corrections to the Wbb process 

pp^ W{-^ eu)g*(^ bb) , (2) 

which is the principal background for the associated Higgs production ([!]) at the Tevatron. 

Other non-negligible backgrounds, provided by the production of WZ, tb, ti and other 

processes ||, will not be considered here. The calculation is performed in the limit 

where 6-quark is massless, and our results indicate that the next-to-leading order (NLO) 

corrections to the Wbb background are significant. 

1 We note here that in the minimal supersymmetric extension of the standard model (MSSM) the Higgs 
boson mass is constrained to be less than about 130 GeV ||, and that for some range of parameters the 
neutral Higgs boson of the MSSM behaves like the Higgs boson of the standard model j|. 
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In Section |2| we briefly describe the calculation method based on the subtraction 
procedure |J, as formulated in J7J. Results for the NLO corrections to the lowest order 
processes (|l|) and (@) at the Tevatron are presented in Section |], while conclusions are 
given in Section [|. 

2 Calculation method 

In order to evaluate the strong radiative corrections to processes ([]]) and (|2]) we have to 
consider Feynman diagrams describing real radiation, and also the ones involving virtual 
corrections to the tree level graphs. 

The corrections due to real radiation are dealt with using the general subtraction 
algorithm formulated by Catani and Seymour @, which is based on the fact that the 
singular parts of the QCD matrix elements for real emission can be singled out in a process- 
independent manner. By exploiting this observation one can construct a set of counter- 
terms that cancel all non-integrable singularities appearing in real matrix elements. The 
NLO phase space integration can then be performed numerically in four dimensions. 

The counter-terms that were subtracted from the real matrix elements have to be 
added back and integrated analytically in n dimensions over the phase space of the extra 
emitted parton, leading to poles in e = (n — 4)/2. After combining those poles with the 
ones coming from the virtual graphs all divergences cancel, so that one can safely perform 
the limit e — > and carry out the remaining phase space integration numerically. 

For the signal process (|l|) we consider only the effects of the initial state gluon emission. 
The final state radiation can be taken into account in the total rate by using the radiatively 
corrected branching ratio for H — > bb ||. The virtual corrections to ([!]) are of the Drell- 
Yan type and are well known |§. They are expressible as a multiple of the lowest order 
matrix element squared. 

For the Wbb background process (0) we consider real radiation from both initial and 
final state quarks. The virtual corrections to the tree level graphs can be obtained by 
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crossing the one loop helicity amplitudes for the process e + e — > qqQQ JTO] . 



Note that our final results are presented in the MS renormalization and factorization 
scheme. However, in intermediate steps for the Wbb process we used the four dimensional 
helicity scheme of |L0] . 



For the sake of simplicity, we performed the calculation in the limit where 6-quark is 
considered massless, and with CKM matrix elements V u b and V c b set to zero. In the case 
of Wbb background, the latter approximation eliminates the need to take into account 
scattering processes involving 6-quarks in the initial state. Given that, for example, 
\V c b/V u d\ 2 ~ 0.002, the effects of setting V c b = V u b = are small. Corrections for the finite 
6-quark mass are expected to be of order 4m^/M^, or about 1% for M b i ~ 100 GeV. 

In order to ensure that we have a hard sub-process, we have also introduced a set of 
basic cuts, 

(Pb + Pb? > 4Q 2 , 

Pi > Q, (3) 

Pi > Q, 

where Q is a scale of the same order as the 6-quark mass. The first constraint imposes 
the correct physical threshold even though we have set the b quark mass to zero. The 
constraint on the of the b and b quarks obviates the need for factorization subtractions 
involving the lowest order process qb — > Wq'b. In general, more stringent cuts on all three 
quantities will be required for comparison with experimental data. 



3 Results 

All results given in this paper report on the rate obtained for W + production in pp 
collisions at v^S" = 2 TeV. To include the contributions from W~ production and the 
contributions from the decay into muons, one should multiply our results by a factor of 
four. Note however that we assume perfect efficiency e& for detection of 6-jets. Achievable 
values of this efficiency would decrease our results by a factor of e\ ~ 0.2 ||TT| . We used 
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the MRSR2 parton distribution functions with as{Mz) = 0.120 [^2] , while the scale Q 
from (^) was set equal to 4.62 GeV. 

In Figure [l] we first show the scale dependence of the LO and NLO cross sections 
for the signal and background processes in a mass window 84 < M bb < 117 GeV, which 
is appropriate Q for a 100 GeV Higgs boson. No other cuts apart from (|3|) have been 
applied. At a natural scale of \i = 100 GeV we find that the K-factor is about 1.2 for the 
signal, but 1.5 for the background. 

As already mentioned, more strict cuts than those of @ have to be applied for com- 
parison with experiment. In addition to the cuts on rapidity and transverse momentum, 



(4) 



we also impose isolation cuts, 

R bb , R eb , R eb > 0.7 , (5) 

as well as a cut on the scattering angle of the bb system |13j (the Higgs scattering angle) 
in the Collins-Soper frame 

|cos%|<0.8. (6) 



\Vb\, \yi\ 


< 


2, 


\Ve\ 


< 


2.5, 


\Pbl \Ph\ 


> 


15 GeV , 


\Pel \Pl\ 


> 


20 GeV , 



Note that imposing the cut on cos 9 bb requires knowledge of the longitudinal component of 
a neutrino momentum. By assuming that W boson is on shell, and using p e and which 
are actually measured, this component can be reconstructed up to a two-fold ambiguity 
for a solution of a quadratic equation. Due to the asymmetry of the neutrino rapidity 
distribution, by choosing the larger (smaller) solution for p z v in the case of W + {W~), 
one can improve the probability of finding the correct W momentum. Following this 
prescription, in our LO (NLO) Monte Carlo studies with the above cuts we have observed 
efficiency of about 77% (68%) and 70% (61%) for processes ([]]) and (|j), respectively.^ 
2 The actual efficiency of accepting an event based on the calculated value of cos 9 b i from the recon- 
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Figures |2|, [| and f| show our results obtained after cuts (J|)-(§D have been imposed. 
In Figure ^ we illustrate the bb mass dependence of the Wbb process in LO and NLO. 
The shape of this curve could be important in extrapolating the Wbb background from 
observed events with lower bb invariant mass. Note however that the we have not yet 
clustered the 6-partons with emitted radiation to form fe-jets, although our numerical 
programs are flexible enough to allow that. This procedure could alter the shape of the 
bb mass spectrum. 

Figures |3] and [| show the scale dependence of the LO and NLO cross sections for 
84 < M bh < 117 GeV (M H = 100 GeV) and 102 < M bh < 141 GeV (M H = 120 GeV). 
At natural scales we estimate a signal to Wbb background ratio of 0.43 and 0.30 at 
Mn — 100 GeV and Mu = 120 GeV, respectively. Before definitive conclusions can be 
drawn on the Tevatron search for the SM Higgs boson other backgrounds need to be 
included. 



4 Conclusions 

We have presented first results from a calculation of the radiative corrections to the pro- 
duction of bb in association with a W. The corrections are observed to be large and 
positive at natural scales. The full implications of this for the search for the standard 
model Higgs boson at the Tevatron will be discussed elsewhere |T^]. Note that sizable 
corrections have also been found in other reactions where the same virtual matrix ele- 



ments have been used in a crossed channel [|16|, [LT|| . As yet we have no clear analytic 
understanding of why this should be so. 

Our calculation can be considered the first step in the calculation of W + 2 jet cross 
section. In addition, extension of our programs to other two boson processes should be 
straightforward. 



structed W momentum is higher. With cuts (H)-® we found that about 93% (87%) and 87% (79%) of 
the LO (NLO) events for processes ([j]) and (||) were correctly accepted or rejected. 
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Figure 1: Cross section scale dependence for the signal (lower curves) and Wbb 
background (upper curves). These results were obtained for 84 < M b i < 117 GeV 
(Mh = 100 GeV). Apart from (|), no other cuts have been applied. 
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Figure 2: The LO and NLO M b i dependence of the Wbb background process at a scale 
of fi = 25 GeV. The results shown were obtained after imposing cuts given in @-(|6]). 
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Figure 3: Cross section scale dependence for the signal (lower curves) and Wbb 
background (upper curves). These results were obtained for 84 < M b i < 117 GeV 
(Mh = 100 GeV), and after imposing cuts given in 
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Figure 4: Cross section scale dependence for the signal (lower curves) and Wbb 
background (upper curves). These results were obtained for 102 < M h i < 141 GeV 
(Mh = 120 GeV), and after imposing cuts given in ©-(0). 
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